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A simple model is proposed to study the sound wave in hot dense matter created in heavy ion 
collisions by jet. The preliminary data of jet shape analysis of PHENIX Collaboration for all 
centralities and two directions is well explained in this model. Then the wavelength of the sound 
wave, the natural frequency of the hot dense matter and the speed of sound wave are estimated 
from the fit. 
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I. INTRODUCTION 

The energy loss of hard partons induced at the early 
stage of heavy ion collision has been regarded as a 
promising tool to study the properties of deconfined 
matter that created in the untrarelativistic heavy ion 
reactions 0 0 0 0, IE IE 0 • While the energy and mo¬ 
mentum that lost by a hard parton propagating through 
hot and dense matter would obviously be redistributed 
in the created medium, recent resultsQ of the mea¬ 
surements of two-particle correlations, which involve on 
hard trigger particle, have shown us a surprising fea¬ 
tures that double peaks are measured on the away side 
for all non-peripheral collisions for Au+Au collisions. 
This distinguished feature is quite different from that 
for Au+Au peripheral collisions and p-p or d-Au colli¬ 
sions. This surprising feature attracts interests of many 
physicists and a lot of models have been proposed to 
explain this phenomenon, e.g., possible colored sound 
mode0 0 EE O E2 proposed by J. Ruppert and B. 
Muller et al(see[U3 for detailed review by J. Ruppert). 

Among these attempts, the idea of conical flow 'Tij . 
similar to sonic booms induced by supersonic bodies, is 
regarded as a promising explanation. In fact, this idea 
has been proposed before the double peaks are observed 
definitely and it predicts this distinguished phenomenon. 
The most excellent property of this model is introducing 
a collective exciting mode to the created medium after 
heavy ion collision. In conical flow model, in order to 
separate the conical flow and the elliptic flow, it is con¬ 
venient to focus on the most central collisions, where the 
elliptic flow is as small as possible. The formalism used 
in this model is based on linearized hydrodynamics,which 
holds true in the case the total energy-momentum den¬ 
sity deposited by the jets is a small perturbation com¬ 
pared to the total energy of the medium when describing 
medium far enough from the jet head.The assumptions 
that the perturbed medium is homogeneous and at rest 
are also required. Although comparing with the angular 
correlation pattern of hard hadrons shows good agree¬ 
ment with theory and experiment,the shortcoming of the 
conical flow model is obvious. 


To study dependence of the angular correlation pat¬ 
tern to centrality i.e. the impact parameter,this letter 
derives the merit of conical flow model and uses a more 
microscopic and basic formalism of this effect. The col¬ 
lective exciting mode means the perturbation induced by 
the jet is not separately. There is correlation between 
contiguous perturbations. In our model, we introduce 
the intrinsic sound wave of relativistic hot dense mat¬ 
ter which will be discussed in detail in the following to 
express the correlation between perturbations. 

The letter is organized as follows: In Section EH we 
explain basic framework and assumptions in our model. 
In Section mu this result of the model is discussed and 
compared with experiment data. A summary is made in 
Section IIVI 


II. FRAMEWORK AND ASSUMPTIONS 

For the shortage of knowledge of interaction between 
parton and gluon, it is difficult to know the behavior 
of hot and dense matter created in heavy ion collision. 
However, to study the energy losses of fast quenched par- 
ton, one treats two different types of cases separately : 

(i) the radiative losses, producing relativistic gluons; and 

(ii) the scattering/ionization losses, which deposit energy 
and momentum directly into the medium. Both of them 
have relations to the property of the hot dense matter 
created in heavy ion collision. 

To overcome this shortage, one can understand this 
energy transportation on the point of view of a wave. 
Generally, a wave can travel through whatever matter 
which satisfies dP/dV < 0 with P and V being pressure 
and volume of the matter respectively,like lattice wave in 
solid, sound wave in liquid and gas,carrying energy. 

It is reasonable to deem dP/dV < 0 is true for hot 
dense matter produced in heavy ion collision since it is 
wrong only for undetected dark matter which is impor¬ 
tant in cosmology. If a wave travelling the matter freely, 
the velocity and frequency of it are intrinsic properties 
of the matter. The hot dense matter created in heavy 
ion collision can be regarded as a kind of liquid or gas, 
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so we call the wave travelling through it as sound wave. 
And the perturbation induced by fast parton will prop¬ 
agate freely until it reaches the boundary of the source. 
So the frequency the sound wave should be the natural 
frequency of the hot dense matter. 

As a standard wave equation always successes in de¬ 
scribing wave phenomena in having equilibrated system, 
we assume it to be able to describe the sound wave in 
hot dense matter for a admittedly thermalization can be 
reached after heavy ion collision. The speed of sound c s 
is the unique parameter of the wave equation and one pa¬ 
rameter in our model which should be invariant in every 
case of centricity. Then sound wave travelling the hot 
dense matter obeys a wave equation reading 

^2 ~ c 2 s Au = 0 W 

where c s is the speed of sound of the medium produced 
in heavy ion collision. We note that equation m is only 
a approximation of true one without considering anhar- 
monic terms, which is reasonable where the anharmonic 
terms inducing dissipative effect in the medium is negligi¬ 
ble. This approximation is also necessary for Mach-Cone 
description and its reasonability should be tested in fur¬ 
ther study of non-perturbation QCD. However,equation 
m at least holds true where is away from the trajectory 
of quenched fast parton for a distance of order of sound 
attenuation length ,T = (4/3)r]/(e+p), with ry being the 
shear viscosity. 

On the other hand ,as the quenched fast parton is con¬ 
stantly emitting gluons, which emit new ones etc., the 
whole shower is a complicated nonlinear phenomenon 
Eii |. The multiplicity of this shower grows nonlinearly 
with time. Eventually, the core may become a macro¬ 
scopic co-moving body. All perturbations originating 
from this co-moving system produced by the quenched 
fast parton requires the consideration of Doppler effect. 
The velocity of the co-moving system c c is another pa¬ 
rameter of our model which should be determined in ex¬ 
periment fit. However, it can also be discussed quali¬ 
tatively. c c should increase with the increment of the 
quenched parton and decrease with the increment of 
the length of the hot dense matter through which the 
quenched parton travels. We note that the source is ex¬ 
panding instead of static,so the whole medium through 
which sound wave travels has a radial or elliptical ve¬ 
locity which should affect the energy distribution of the 
final-state particles. This effect is approximately negli¬ 
gible in center-center collision because a radial velocity 
field will not change the spatial angle distribution of en¬ 
ergy. For a simple case, we adopt a radial velocity field 
in our discussion. 

Now we depict our image of understanding the interfer¬ 
ence of perturbation induced by the quenched fast parton 
in Fig.l, where perturbations originate from line section 
AB with phase 2nlc s /(VX) where l is the distance be¬ 
tween B and point located in line section AB, V is the 
velocity of the quenched fast parton ,always taking c for 



FIG. 1: diagrammatic sketch of perturbation interference 
induced by quenched fast parton. perturbations originates 
from A and B has a phase difference 2nL-g(j/\ where A is 
the wavelength of sound wave traveling through the medium 
created in heavy ion collision 

it, and A is the wavelength of the sound wave. The per¬ 
turbations interfere with phase and produce the inter¬ 
ference pattern ,for example, in D and E. This pattern 
represents the distribution of energy of perturbations and 
should be equivalent in shape to two-particle correlation 
about jet which is measured in experiment. 

The momentum distribution of particles induced by 
the quenched fast parton is not treated in this model 
since it has relations to durative scattering rather than 
collective excitation mode. For this reason, we can only 
adjust the relative intensity between interference pattern 
of trigger and its counterpart according to experiment. 

Conclusively, there are three parameters in our model 
totally, which are the speed of sound c s . the velocity of 
co-moving system c c and the ratio of the size of the source 
to the wavelength of the sound wave Lj^/ A. Because the 
system in which we are interested evolutes with a finite 
lifetime, these parameters should be regarded as time- 
weighted averages respectively. 


III. RESULT 

In experiment, one can measure the azimuthal correla¬ 
tion ,around a high momentum trigger, of particles within 
a special momentum window , which should be repro¬ 
duced by the diffraction pattern of the carrying energy 
sound wave in our model. 

The perturbations induced by a same parton should be 



The parameters in the calculation is listed below 
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coherent by phase, while these induced by different par- 
tons should be superposed by intensity. Thus the diffrac¬ 
tion pattern for one parton reads 


1(6) = A 


sin[ 7 rs(cos 0 — c s )/A] 2 
7 rs(cos 9 — c s )/A 


( 2 ) 
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where 9 is the azimuth angle in the frame of center of 
mass and s is the length through which the fast parton 
passes and A is the wavelength of the sound wave. For 
showing the result in laboratory frame, one can substi¬ 
tute 9' for 9 with, 


9’ = 


f tan 1 | 

( p sin 6 \ 

1 7(pcos 6—(3E) J 

) tan -1 | 

{ p sin 6 \ 

^7(pcos 6—(3E) J 


Sign(9)n 1 


pcos9 > /3E, 

otherwise. 


(3) 


where /3 is c c and 7 = (1 — (3 2 )( 2 ). p and E in m 
means the momentum and energy of observed particle in 
laboratory frame. The analytic form of 1(6') is difficult 
to write out, though it looks very easy. 

The speed of sound c s should be a constant for a fixed 
temperature and the speed of co-moving system c c should 
increase obviously with the momentum of trigger par¬ 
ticle. the length of trigger passing s should be deter¬ 
mined by the centricity of the collision, i.e. the multi¬ 
plicity of the collision. Fig.l shows the comparison be¬ 
tween prelimi nary data of jet shape analysis of PHENIX 
Collaboration[l5| and our calculation. Different marks 
denote experiment results in different centralities, while 
different styles of lines denote our calculations using dif¬ 
ferent parameters. The momentum of particle used in 
the calculation, 2.5 GeV , is the average of the momen¬ 
tum range of experimentally measured particles i.e. 2-3 
Gev. The used speed of sound , much smaller than that 
used in Mach cone calculation in Hi where c s = 0.33c, 
is 0.129c for the better fit, which is a constant in calcu¬ 
lation. In this case the peak of induced particle should 
be at angles about 1.7 radian relative to the trigger. The 
normalization for all centralities has been performed. 



The | means the forward jet i.e. the jet caused by the 
trigger particle while J, means the backward jet which is 
caused by the counterpart of the trigger. 

It is easy to understand the variation of the parame¬ 
ters. The centrality of the collision , relative to the source 
size, is reflected by the s/ A. As the increment of the 
source size, the co-moving velocity decreases.This reflects 
the energy loss of particle when traveling the source. 

The intensity ratios between forward jet and backward 
jet are 0.60, 0.51 and 0.67 for 0-10%, 30-40%, 60-92% 
respectively, which come from data fitting and will not 
be discussed in our model as is above mentioned. 

Now the picture of sound wave in hot dense matter is to 
some extend obvious. Let’s take double HBT radius 
, 5 /to, as the source size for center-center collision, the 
wavelength of sound wave in hot dense matter created in 
heavy ion collision is estimated as A = (2 x 5/m)/(0.6 + 
0.15) = 13.3 fm. The natural frequency of the hot dense 
matter is estimated as c s /A = 2.91 x 10 21 Hz. These 
quantities will be constants at fixed temperature. 

But Making conclusion that we have catch signals 
about the intrinsic eigenvibration of hot dense matter 
is premature at present. After all, there must be other 
models being able to explain the same data. To test the 
existence of the sound wave, we can examine other be¬ 
haviors of this model. 

The velocity of co-moving system c c will increase as the 
increment of the momentum of the trigger. For higher 
momentum of trigger particle,the distance between two 
peaks in backward jet in center-center collisions will de¬ 
crease since the Lorentz transformation in ®), which can 
be tested in LHC where a lot of jets with high momentum 
trigger will be produced. 

We can also make a more fine bins of centrality to draw 
the jet induced particle azimuth correlation again, which 
will give a more crashing test of this model. 

Essentially, the properties of the sound wave are deter¬ 
mined by the interaction between particles in hot dense 
matter. One can conceive that it is quite different for 
confined and deconfined scenarios. So it is meaningful to 
compare different scenarios to experiment, which will be 
a interesting tool to research QGP. 


IV. CONCLUSION 


FIG. 2: the comparison between experiment and theory, The 
dash line corresponds to center-center collisions, where the 
size of source is relative large. The thin corresponds to pe¬ 
ripheral collisions. The particle momentum is 2.5 Gev and 
c s = 0.129 


We try to study the sound wave in hot dense matter 
created in heavy ion collisions in this letter by jet. The 
distribution of final-state particle is described by the in¬ 
terference of the sound wave induce by the quenched fast 
parton. The productions of a co-moving system around 
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the quenched fast parton is introduced. The preliminary 
data of jet shape analysis of PHENIX Collaboration for 
all centralities and two directions is well explained in our 
model . From the parameter of the result, we estimate 
the wavelength of the sound wave is 13.3 fm, the natural 
frequency of the hot dense matter is 2.91 x 10 21 Hz and 
the speed of sound wave is 0.129c. 

We predict that the distance between two peaks in 
backward jet in center-center collisions will decrease since 
the Lorentz transformation for higher momentum of trig¬ 
ger particle, which can be tested in LHC where a lot of 
jets with high momentum trigger will be produced in the 
future. Our model cares only collective property of the 


hot dense matter. Whether the sound wave can become 
a tool to test the creation of quark-gluon plasma depends 
on the difference between calculations of the sound wave 
in quark-gluon plasma scenario and hadron-gas scenario. 
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